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Abstract 

DB boundary conditions have been simulated in terms of perfect electric conductor 
(PEC) boundary for transverse electric modes and perfect magnetic conductor (PMC) 
for transverse magnetic modes. Electric and magnetic fields inside the dielectric region 
of a parallel plate DB boundary waveguide have been derived and fractional curl 
operator has been utilized to study the fractional parallel plate DB waveguide. Field 
behavior as well as transverse impedance of the guide walls have been studied with 
respect to the fractional parameter describing the order of fractional curl operator. 
The results are also compared with the corresponding results for fractional waveguides 
with PEC walls. 
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1. Introduction 

Fractional calculus is a branch of mathematical analysis that studies the possi- 
bility of taking real number powers or complex number powers of the differentiation 
and integration operators [1]. It fills in the gaps of classical calculus, so in this re- 
spect traditional calculus may be taken as a special case of the fractional calculus. 
Using fractional calculus, scientists and engineers have been interested in exploring 
the potential utilities and possible physical implications of mathematical machinery 
of the subject of fractional calculus, that is, fractional derivatives and fractional in- 
tegrals [2]. Fractional derivatives and integrals are mathematical operators involving 
differentiation and integration of arbitrary (non-integer) real or complex orders such 
that d'^ f (x) / dx'^ and d~'^ f (x) / dx~'^ respectively, where m and n can be taken as 
non- integer real or even complex number. It has been demonstrated that these math- 
ematical operators are interesting and useful tools in various disciplines of science and 
engineering [3-5]. 

Electromagnetic theory has a driving role in the modern world of science and 
engineering. Maxwell equations encapsulate all this revolutionary discipline, whose 
solutions have a great importance in current research and development activities. Fif- 
teen years ago, Engheta particularly focused on finding out what possible applications 
and/or physical role, the mathematical operators of fractional calculus can have in elec- 
tromagnetic theory [6-11]. He applied the concept of fractional derivatives/integrals 
to certain electromagnetic problems, and obtained interesting results and ideas show- 
ing that these mathematical operators are interesting and useful mathematical tools in 
electromagnetic theory. Some of these ideas include the mathematical link between the 
electrostatic image methods for the conducting sphere and the dielectric sphere [6], 
fractional solutions for the scalar Helmholtz equation [7,10], electrostatic fractional 
image methods for perfectly conducting wedges and cones [8], and the novel concept 
of fractional multipoles in electromagnetism [9] . 
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Other researchers belonging to the area of electromagnetics have also contributed 
to this topic. Tarasov proved that the electromagnetic fields in dielectric media, whose 
susceptibility follows a fractional power-law dependence in a wide frequency range, can 
be described by difi^erential equations with time derivatives of noninteger order [12,13]. 
Fractional dimensional space represents an effective physical description of confinement 
in low- dimensional systems. The concept of fractional calculus to obtain the solution 
of electrostatic problem in fractional dimensional space, for the fractional order mul- 
tipoles, is utilized by Muslih and Baleanu [14]. They also introduced the form of 
fractional scalar potential by using the solutions of Laplace's equation in fractional di- 
mensional space [15]. They derived potential of charge distribution in fractional space 
using Gegenbauer polynomials. According to Zubair et al., solutions of Helmholtz 
equation in fractional space can describe the complex phenomenon of wave propagation 
in fractal media. With this view, they established a generalized Helmholtz equation 
for wave propagation in fractional space and found its analytical solution [16-21]. 

Fractionalization of ordinary derivative and integral operators motivated the re- 
searchers to find out the possible fractionalization of other operators and their use in 
electromagnetics. Engheta fractionalized the kernel of an integral transform and stud- 
ied the paradigm of intermediate zones in electromagnetism [22,23]. Fractionalization 
of the curl operator, a well known operator in vector calculus, was also introduced 
by him [24]. He used the fractional curl operator to find the new solutions, which 
may be regarded as intermediate step between the two given solutions, to the Maxwell 
equations. 

For the sake of completeness, it is decided first to reproduce the meanings of the 
following concepts given by Engheta [24] : what is meant by fractionalization of a linear 
operator? How to fractionalize a linear operator? That is, what is the mathematical 
recipe to fractionalize a linear operator? Fractionalization of the curl operator using 
this recipe and utilization of the fractional curl operator to find intermediate/fractional 



-4- 

solutions to the Maxwell equations for ordinary medium has been discussed [24] and 
just results of their work are presented here. Time harmonic dependency exp[—iojt) 
has been considered throughout the paper. 

1.1. Fractional linear operator in electromagnetics 

1.1.1. Conditions required to verify fractionalization 

Mathematical fractionalization of any problem requires two canonical solutions of 
the problem under consideration and an operator that can transform one canonical 
solution into the other. Fractionalization of the connecting operator can reveal inter- 
mediate solutions between the two canonical solutions. The conditions and recipe for 
fractionalization of a linear operator L are reproduced below [24]. It must be men- 
tioned that this recipe has also been used to fractionalize the Fourier transform [25]. 
The new fractionalized operator L'* with fractional parameter a, under certain condi- 
tions, can be used to obtain the intermediate cases between the canonical case 1 and 
case 2. 

A linear operator L may be a fractional operator (i.e., L") that provides the inter- 
mediate solutions to the original problems, if it satisfies the following properties [11,24]. 

I. For Q! = 1, the fractional operator L" should become the original operator L, 
which provides us with case 2 when it is applied to case 1. 

II. For a = 0, the operator If should become the identity operator I and thus 
the case 1 can be mapped into itself. 

III. For any two values ai and a2 of fractional parameter, L'* have the additive 
property in a, i.e.. 



- 5 - 

IV. The operator should commute with the operator involved in the mathe- 
matical description of the original problem. 

1.1.2. Recipe for fractionalization 

It is assumed that, present discussion is about a class of linear operators (or 
mappings) where the domain and range of any linear operator of this class are similar 
to each other and have the same dimensions. That is, L'* : — >• where C"^ 
is a n dimensional vector space over the field of complex numbers. Once a linear 
operator such as L is given, the recipe for constructing the fractional operator can 
be described as follows [24]. 

1. One finds the eigenvectors and eigenvalues of the operator L in the space C"^ so 
that L.Ato = CLm-^m where and am for m = l,2,3....n, are the eigenvectors 
and eigenvalues of the operator L in space C"' respectively. 

2. Provided A^s form a complete orthogonal basis in the space C"', any vector in this 
space can be expressed in terms of linear combination of A^. Thus an arbitrary 
vector G in space C" can be written as 



where gmS are co-efficients of expansion of G in terms of A^s. 

3. Having obtained the eigenvectors and eigenvalues of the operator L, the frac- 
tional operator can be seen to have the same eigenvectors A^s but with the 
eigenvalues as (am)'^, i.e.. 



When this fractional operator operates on an arbitrary vector G in the space 
C"^, one gets 




(1.2) 



m 



(1.3) 



L".G = L"V" OmA 




(1.4) 
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The above equation essentially defines the fractional operator L'* fi:om the knowl- 
edge of operator L and its eigenvectors and eigenvalues. In the next section, above 
recipe has been applied to fractionalize the curl operator. 

1.1.3. Fractional curl operator and Maxwell equations 

Consider a three-dimensional vector field F as a function of three cartesian space 
coordinates {x,y,z). Curl of this vector can be written as 

curlF - (^^^ (^^"^ _ V (^^'^ - ^^^^ z (1 5) 

\ dy dz \ dz dx J ^ \ dx dy J 

where Fx, Fy, F^ are the cartesian components of vector F and x, y, z are the unit 

vectors in the space domain. Assuming that spatial Fourier transforms of both the 

vector functions (F and cur IF) exist, the Fourier transform of these two vectors can 

be written as 

J^k {F(a;, 2/, z)} = F(/c^, ky, k^) 

/OO I'OO POO 
/ / F{x,y, z) exp{—ikxX — ikyy — ikgz)dx dy dz{lS) 
-OO J —QO J —OO 

/OO poo poo 
I I cmlF{x, y, z) exp{—ikxX — ikyy — ikgz)dx dy dz 
-OO J —OO J —OO 

= i'kxF{kx,ky,k^) (1.7) 

where a tilde over the vector F denotes the Fourier transform of vector F. Hence in 
the k-domain {kx, ky, kz), the curl operator can be written as a cross product of vector 
zk with the vector F. In order to fractionalize the curl operator, the cross product 
operator (zkx) in the k-domain serves as prerequisite. Thus fractionalization of curl 
operator is equivalent to the fractionalization of this cross product operator. With the 
recipe described in the previous section, fractionalization of the cross product operator 
as {ikx)°' can be obtained in the k-domain. 

Engheta utilized the fractional curl operator to fractionalize the principle of du- 
ality in electromagnetics [24] . Fractionalization of the principle of duality yields new 



solutions to the Maxwell equations, which may be regarded as intermediate step be- 
tween the original solution and dual to the original solution and has been termed as 
fractional dual solutions. In an isotropic, homogeneous, and source free medium de- 
scribed by wavenumber k and impedance rj, new set of solutions to the source-free 
Maxwell equations may be obtained using the following relations [24] 



Efd = 



7Ti^(ikx)"E 
-T^(ikx)"?7H 



(1.8a) 
(1.85) 



where fd stands for the fractional dual. Inverse Fourier transforming these back into 
the {x, y, 2)-domain, the new set of solutions are obtained as 



Efd = 



1 



{iky 

1 

Jik)' 



-curl^E 



(1.9a) 
(1.96) 



From Eqs. (1.9), it can be seen that for a = 0, (Efd, r/Hfd) gives the original solutions 
whereas (Efd, r/Hfd) gives dual to the original solution to the Maxwell equations for 
a = 1. Therefore for all values of a between zero and unity, (Efd,?7Hfd) provides 
the new set of solutions which can effectively be regarded as intermediate solutions 
between the the original solution and dual to original solution. These solutions are 
also called the fractional dual fields as expressed with the subscript fd. 

1.4. Previous contributions 



As it has already been stated that concept of the fractional curl operator and its 
utilization in electromagnetics was given by Engheta [24]. Naqvi and Rizvi extended 
Engheta's work on fractional curl operator by determining sources corresponding the 
fractional dual solutions to the Maxwell equations. Results of their valuable work show 
that surface impedance of the planar reflector, which is intermediate step between the 
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PEC and PMC, is anisotropic in nature [26]. Naqvi et al., further extended work on 
this topic by finding fractional dual solutions to the Maxwell equations for reciprocal, 
homogenous, and lossless chiral medium [27]. Lakhtakia pointed out that any frac- 
tional operator that commutes with curl operator may yield fractional solutions [23]. 
Naqvi and Abbas studied the role of complex and higher order fractional curl operators 
in electromagnetic wave propagation [28] . They also studied the fractional dual solu- 
tions in double negative (DNG) medium [29]. Veliev further extended the work on the 
fractional curl operator by finding the refiection coefficients and surface impedance 
corresponding to fractional dual planar surfaces with planar impedance surface as 
original problem [30]. The work on this topic entered into new era when concept of 
fractional transmission lines, fractional waveguides, and fractional resonator in elec- 
tromagnetics were introduced [31-40] and nature of the modes supported by fractional 
dual waveguides and impedance of the walls were addressed. Modelling of transmission 
of electromagnetic plane wave through a chiral slab using fractional curl operator and 
fractional dual solutions in bi- isotropic medium are also available [41,42] 

After the introduction of nihility concept by Lakhtakia [43], Tretyakov incorpo- 
rated the nihility conditions to chiral medium and proposed another metamaterial 
termed as chiral nihility metamaterial [44,45]. Chiral nihility is a metamaterial with 
following properties of constitutive parameters at certain frequency [45]. 



Thus the resulting constitutive relations for isotropic chiral nihility metamaterial re- 
duces to 



e ^ 0, 



(1.15a) 



B = - 



(1.156) 
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Tellegen nihility [46] states that 

e^O, fJ'-^O, K^O, Xt^O 

and corresponding expressions for constitutive relations for Tellegen nihility metama- 
terial are 



Study of nihility /chiral nihility metamaterials is a topic of current research by several 
researchers [47-57]. Naqvi contributed many research articles on chiral nihility and 
fractional dual solutions in chiral nihility metamaterial [51-57]. 

1.3. DB boundary conditions 

Before the advent of idea 'DB boundary interface' proposed by Lindell and Si- 
hvola[58,59], all the known interfaces dealt with tangential components of electric and 
magnetic fields. But the DB boundary is analyzed on the basis of normal components 
of flux densities D and B [ 59,60]. Waves polarized transverse electric (TE) and trans- 
verse magnetic (TM) with respect to the normal of the boundary are reflected as from 
perfect electric conductor (PEC) and perfect magnetic conductor (PMC) planes ( i.e., 
DB interface behaves like PEC and PMC) [61]. It is worth mentioning here that, all 
of the previous conditions in electromagnetics are associated to the electromagnetic 
fleld vectors E and H. The boundary conditions for DB interface may be written as 



D = xV^o/^oH 



(1.16a) 



(1.166) 



[59-65] 



n.D = 



(1.17a) 



fi.B = 







(1.176) 



where n is normal vector to the interface. 



Purpose of the current discussion is to extend the previous work on fractional 
waveguides for waveguide with DB walls. 



2. General wave behavior along a guiding structure 

Consider a waveguide consisting of two parallel plates separated by a dielectric 
medium with constitutive parameters e and ix. One plate is located at ?/ = 0, while 
other plate is located at y = b. The plates are assumed to be of infinite extent and 
the direction of propagation is considered as positive z-axis. Electric and magnetic 
fields in the source free dielectric region must satisfy the following homogeneous vector 
Helmholtz equations 

V^E{x,y,z) + k'^E{x,y,z)^0 (2.1a) 
W'^n{x,y,z) + k'^n{x,y,z) = (2.16) 

where V = ^2 + ^2 + ^2 is the Laplacian operator and k = is the wave 

number. Taking z-dependance as exp{iPz), equation (2.1) can be reduced to two 
dimensional vector Helmholtz equation as 

VlyE{x,y) + h^E{x,y) = (2.2a) 
Vlynix,y) + h^liix,y) = (2.26) 

where h'^ = k'^ — /3'^ , /3 is the propagation constant. 

Since propagation is along z-direction and the waveguide dimensions are consid- 
ered infinite in xz-plane. So x-dependence can be ignored in the above equations. 
Under this condition, equation (2.2) becomes ordinary second order differential equa- 
tion as 

+ h'E{y) = (2.3a) 
+ h^u{y) = (2.36) 

ay 

As a general procedure to solve waveguide problems, the Helmholtz equation is solved 
for the axial field components only. The transverse field components may be obtained 
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using the axial components of the fields and Maxwell equations. So scalar Helmholtz 
equations for the axial components can be written as 

"^'^^ +h'E, = Q (2.3c) 



+ h'^H^ = {2.3d) 



dy 

General solution of the above equations is 

Eg = an cos(hy) + bn sm(hy) (2.3e) 

Hg = Cn cos{hy) + dn sm{hy) (2.3/) 

where an,bn,Cn, and dn are constants and can be found from the boundary conditions. 

Using Maxwell curl equations, the transverse components can be expressed in 
terms of longitudinal components {E^, Hg), that is 

^ h"^ dx ' rj dy ) (^•^'^) 
1 f .r,dHg ikdEA 



where 



\ dy T] dx 



ry = Y — is impedence of the medium inside the guide 

In the proceeding part of this paper, parallel plate waveguide with DB boundary- 
walls has been considered and the fractional dual solutions have been determined and 
analyzed. 



3. Parallel plate DB boundary walls waveguide 



A wave of general polarization propagating in positive z-direction through a par- 
allel plate waveguide can be written as a linear sum of the transverse electric(T£'^) 
modes and transverse magnetic (TM^) modes. A DB boundary can be simulated as 
the boundary which behaves like perfect electric conductor (PEC) boundary for (TE^) 
modes and perfect magnetic conductor (PMC) boundary for (TM^ )modes. Therefore 
fields inside a parallel plate DB waveguide may be obtained by linear superposition 
of two canonical solutions which are transverse electric (TE*^) mode solution for PEC 
waveguide and transverse magnetic(TM^) mode solution for PMC waveguide. Both 
the cases have been discussed separately. 

3.1 Fractional dual solutions of canonical cases 

Case 1: Transverse electric (TE^) mode propagation through a PEC 
waveguide 

Let us first consider that (TE^) mode is propagating through a PEC waveguide 
described in section 2. For this mode, axial component of the electric field is zero 
while for the magnetic fields it is given as in equation (2.3f). Using equations (2.4), 
the corresponding transverse components can be written as 



Using boundary conditions for PEC boundary that is E^^z = ^\y=o,b > we get the 
particular solutions as 




(3.1a) 




K = 



(3.1c) 



H, = 



(S.ld) 




(3.2a) 
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TT 1 ^ I r y . / J \ 1 

vHy = ( y ) [-Cnsm(%)j 




(3.26) 


r/iy^ = CnCOs{hy) 




(3.2c) 


Ey^O 




{3.2d) 


H, = 




(3.2e) 


Cn — Cnf] h = —r- 


n = 1,2,3... 





where 



Re-introducing the z-dependance exp{iPz) and writing equations 3.2 in exponential 
form we can write electric and magnetic fields inside the dielectric region as sum of 
two plan waves as 

E = El + E2 (3.3a) 
r/H = r/Ei + r/Hs (3.36) 

where (Ei,Hi) are the electric and magnetic fields associated with one plane wave 
and are given below 

El = f f ^x) exp{-ihy + i^z) (3.4a) 



2 J \h 

ryHi = ( ^ j (^^y + z ) eM-ihy + i^z) (3.46) 

while (E2, T]iiL2) are the electric and magnetic fields associated with the second plane 
wave and are given below 

E2 = (~^'^) ®^P(^^^ + ^'^-^) (^•^") 

r/Hs = (-^y + z) ex^iihy + i^z) (3.56) 
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This situation can be shown as in Figure 1. 



,=& ///////////////////////////////// 




/77777777777777777777777777777777r^ 



Figure 1 Plane wave representation of the fields inside the waveguide 

Once we have written electric and magnetic fields inside the dielectric region in terms 
of two plan waves, recipe for fictionalization [24,32] can be applied to get the fractional 
dual solutions as 



^PECfd ~ I ^ 



-iCo, sm + — j X + -Soc cos + — j y 



r/H 



TE _ ^ 
PECfd — 



h 



i-^Q sm [hy + —jz exp {^(Jz + )\ (3.6a) 



-So, cos (^/iy + — j X - z-Ca sm (^/iy + — j y 



+^Cq cos (^/ly + ^) z exp i (^/3^ + (3.66) 



where 



/Q!7r\ 
= rns 



C„ = COS 



V 2 / 



S(y — sin 



/CtTT 
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Case 2: Transverse magnetic (TM^) mode propagation through a 
PMC waveguide 



Similar to the treatment done in Case 1, using equation (2.3e) and equations (2.4), 
we can write the results for transverse magnetic mode propagating through a PMC 
waveguide as 



^PMCfd — 



-ba COS + — J X - z-Co, sm \hy + —jy 



■ In 



^ttTM _ a I 
7-n-pMCfd — I — 



+ -6aCos [hy + —jz 
iCa sin (hy + x - ^Sa cos (hy + y 



+ -S^s.n[hy+—) 



exp 



i [Pz + 



an 



(3.6a) 



(3.66) 



3.2 Fractional dual parallel plate DB waveguide 

Fractional dual solutions for the DB waveguide can be written by taking linear 
sum of the fractional dual fields of the above two cases as 

^id — -EJpECfd + ^PMCfd 
'7Hfd = r/Hplcfd + 'yHpMCfd 



which give 



Efd = 



-J exp 



i + {"(^n'S'aCy-i-Q; + ^C'nC'a'S'y+Q;)- 



fd 



-jexp 



+ ^{AnCaCy-\.a — iCnSaSy+a)^} (3.7a) 



i \^Pz -\ —J | — (C^/S^Cy-i-Q, — iAjiCaSy-^ce)^t 
~'^{AnSaCy+a + ^C'nC'a'S'j/+a)y 

h 

+ ^{CnCaCy+a + 'i'AnSaSy+a)^} {3.7b) 
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where 




An,Cn are the constant to be determined from initial conditions. The fields given in 
equation (3.7) have been plotted in Figure 2 for different values of a at an observation 
point {hy,l3z) = (7r/4,7r/4). 




Figure 2 Plots of fractional dual TE^ polarized fields at a point {hy,(3z) — (7r/4,7r/4) 
(a) real parts (b) imaginary parts 



From Figure 2, it can be seen that fractional dual fields satisfy the principle of 
duality, i.e., for ex. = 



and for CK = 1 



-E'fdx = -£^x, 
-E'fdy = Ey, 

Eidz = Ez, 



Efdy = VHy: 

Efdz = vHz, 



vHidy = vHy 



vHfdx = —Ex 
VHidy = ~Ey 
vHfdz = —Ez 



4. Results and discussion 



4.1 Behavior of fields inside the fractional parallel plate DB waveguide 

In order to study the behavior of fields inside the fractional parallel plate DB 
waveguide, electric and magnetic field lines are plotted in the yz-plane and are shown 
in Figure 3. 

These plots are for the mode propagating through the guide at an angle 7r/6 so 
that P/k = cos(7r/6), h/k = sin(7r/6). Initial conditions for both the modes are taken 
as same. Solid lines show the electric as well as magnetic field plots for DB waveguide 
while the fields of PEC waveguide are shown by dashed lines as a reference. From the 
figure we see that there is no normal component of the electric as well as magnetic field 
for q; = 0. This is because the plates of the guide behave as perfect electric conductors 
for transverse electric components while they behave as perfect magnetic conductor 
for transverse magnetic modes. For the reference PEC results, there is no tangential 
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component of the electric filed at the guide surface while magnetic field has no normal 
component. As value of the a increase from 0, normal components of both the fields in 
DB guide starts appearing and become maximum at a = 0.5. After this value normal 
components start decreasing and again become zero at a = 1. 




Figure 3 Field lines in yz- plane at different values of a; solid lines are for the 
fractional DB waveguides while dashed lines are for the fractional PEC waveguides 



- 19 - 

It may be noted that electric field distribution in the DB waveguide is same as 
the magnetic field distribution for the limiting values of a while it is different for the 
intermediate values. Further it may be noted that field behavior for the original and 
dual situation is similar. 

4.2 Transverse impedance of walls 

Wave impedance is defined by ratio of the transverse components of the electric 
and magnetic fields as 

^ Efdx k AnSctCy-\-ct + iCnCotSy-\-a. 

-Hfdz h CnCctCyj^a + i^nSctSy.^^ 



^fdzx =-El—-r]T 



a. 
a 



At J/ = 0, these impedances become impedance of the new reflecting boundary called 
the fractional dual boundary. The normalized impedance matrix of the DB boundary 
wall taking — C„ can be written as 

{k _ h _1 
=fd ^ I ^^fdxzxz + -2;fdzxzx > , < a < 1 



where 



■^fdzx — 



CctCa '^SctSfy 



These impedance components have been plotted for whole range of a as in Figure 4. 
It may be noted that ^^fdxz corresponds to the impedance for TE^ modes and -Zfdzx 
corresponds to the impedance for TM^ modes. Since DB boundary behaves as PEC 
for the TE^ modes so it is zero at a = and a = 1 while it is complex for the 
intermediate range of a. Similarly Since DB boundary behaves as PMC for the TM^ 
modes so it is infinitely high at a = and a — 1 while it is finite complex for the 
intermediate range of a. 
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Figure 4 Transverse impedance of walls of fractional DB waveguides vs. fractional 
parameter a 

5. Conclusions 

It may be concluded from the above results that fractional dual solutions of the 
fields inside a parallel plate DB waveguide may be derived using fractional curl opera- 
tor. The waveguide modelled by such fields may be termed as fractional dual parallel 
plate DB waveguide. It has been noted that electric field distribution in the DB waveg- 
uide is same as the magnetic field distribution for the limiting values of a. while it is 
different for the intermediate values. Fractional dual waveguide may predict the situ- 
ation which is an intermediate step of DB boundary waveguide and its dual situation. 
For limiting cases, transverse impedance of DB wall is zero for transverse electric mode 
and infinitely high for transverse magnetic mode while it is non zero complex value for 
the intermediate situations. 
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